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Abstract 


This  thesis  is  an  investigation  of  the  effects  of  composite 
laminate  orientation  on  the  divergence  problem  of  forward  swept  wings. 
The  aerodynamic  and  structural  properties  are  modeled  with  a  doublet 
lattice  technique  and  finite  element  method,  respectively.  The  static 
equilibrium  equation  is  applied  to  combine  the  two  models .  Th-.;  results 

in  an  eigenvalue  problem  with  the  divergence  dynamic  pressure  as  the 
lowest  eigenvalue. 

A  computer  program  is  written  to  formulate  and  solve  the 
eigenvalue  problem.  It  is  applied  to  two  cases  of  example  wings. 

For  the  first  case,  100%  of  the  fibers  are  varied  from  0°  to  50° 
ahead  of  the  structural  sweep  axis  for  five  different  sweep  angles. 

For  the  second  case,  only  70%  of  the  fibers  are  varied  while  the 
remaining  fibers  'are  oriented  at  90°  and  +  45°  to  the  structural  sweep 
axis . 

The  results  Indicate  that  the  divergence  problem  can  be 
effectively  controlled  and  possibly  eliminated  by  properly  orienting 
the  composite  laminates.  The  divergence  speed  was  found  to  increase 
and  reach  a  maximum  as  the  principal  fiber  direction  is  swept  forward. 
Maximum  values  occur  for  fibers  oriented  between  15°  and  20°  ahead 
of  the  structural  sweep  axis  for  wings  swept  forward  more  than  25°. 
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A  STUDY  OF  THE  USE  OF  LAMINATED  COMPOSITE 


WING  STRUCTURES  FOR  CONTROLLING  THE 
DIVERGENCE  SPEED  OF  FORWARD  SWEPT  WINGS 

I .  Introduction 

Background 

The  use  of  forward  swept  wings  in  aircraft  design  Is  practically 
nonexistant.  This  can  be  attributed  to  the  aeroelastlc  phenomenon 
known  as  divergence.  The  divergence  speed  of  a  wing  can  be  defined 
in  the  following  manner: 

When  a  wing  reaches  its  divergence  speed, 
the  increase  of  aerodynamic  forces  due  to  an 
arbitrary  increase  of  the  twiBt  angle,  is  exactly 
equal  to  the  increase  of  elastic  restoring  forces. 

When  this  speed  is  exceeded,  the  increase  of 
elastic  restoring  forces  required  to  counteract 
the  increased  aerodynamic  forces,  exceeds  the 
elastic  limits  of  the  structure.  This  instability 
can  result  in  drastic  deformations  and  the 
destruction  of  the  wing. 

The  excess  structural  weight  of  conventional  metallic  structures 
required  to  counteract  divergence  has  been  so  great  as  to  prohibit 
the  use  of  forward  swept  wings.  Since  the  divergence  problem  is 
practically  nonexistant  in  conventional  swept  back  designs,  little 
has  been  learned  on  how  to  control  divergence. 
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However,  swept  forward  wings  are  aerodynamically  superior  to  swept 
back  designs.  Several  important  advantages  are  discussed  by  both 
Krone  (Ref  6)  and  Weisshaar  (Ref  7).  The  reports  published  by  these 
two  men  Indicate  that  the  divergence  problem  can  be  effectively 
eliminated  by  the  proper  use  of  advanced  composite  materials.  In 
his  report,  Krone  demonstrated  that  the  use  of  composite  materials 
for  wing  surfaces  could  control  divergence  without  adding  significant 
amounts  of  structural  weight.  Weisshaar  combined  aerodynamic  strip 
theory  with  a  simple  box  structure  to  predict  that  composite  fiber 
orientations  of  10°  to  15°  ahead  of  the  structural  aweep  axis  are  the 
most  effective  in  controlling  divergence. 

It  is  the  purpose  of  this  paper  to  continue  the  study  of 
composite  materials  and  their  affect  on  the  divergence  speed  of 
forward  swept  wings.  An  attempt  will  be  made  to  verify  the  conclusions 
of  Krone  and  Weisshaar  using  methods  different  from  theirs. 


Scope 

An  eigenvalue  problem  was  developed  to  study  the  affects  of 
composite  fiber  orientation  on  the  divergence  speed  of  forward  swept 
wings  A  doublet  lattice  method  was  used  to  model  the  aerodynamic 
forces.  It  represents  the  continuous  lift  distribution  as  a 
discrete  set  of  forces  dependent  upon  the  rigid  and  elastic  orientation 
of  the  winjr..  Similarly,  the  structural  properties  were  modeled  with  a 
finite  ei.em  nt  method.  The  static  equilibrium  equation  was  used  to 
relate  the  structural  deformations  to  the  externally  applied  forces. 
Combination  of  the  two  models  results  in  an  eigenvalue  problem 


2 


with  the  divergence  dynamic  pressure  as  the  lowest  eigenvalue. 


NASTRAN  was  used  to  determine  the  flexibility  matrix  for  the 
structure.  Another  program  combines  NASTRAN' s  flexibility  matrix 
with  an  aerodynamic  matrix  calculated  from  the  doublet  lattice  method. 
It  formulates  the  eigenvalue  problem  and  solves  it  for  the  divergence 
speed.  This  program  was  applied  to  several  example  wings  and  the 
composite  lamina  orientations  were  varied.  Divergence  speeds  were 
calculated  and  compared.  Conclusions  were  made  on  how  the  divergence 
speed  is  affected  by  composite  fiber  orientation. 
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II.  Formulatio  of  Method  for  Calculating  Divergence  Speed 


Introduction 

The  development  of  the  eigenvalue  problem  begins  with  the  basic 
equilibrium  equation: 

F  -  k6  (1) 

This  equation  relates  the  static  aerodynamic  forces  (F)  to  the  wing 
displacements  (6).  The  typical  wing  is  a  very  complex  structure 
consisting  of  many  intricate  structural  members.  Because  of  this  it 
is  all  but  impossible  to  develop  an  accurate  analytic  expression  for 
the  stiffness  coefficient  (k) .  The  resultant  fact  is  that  the 
equilibrium  equation  is  very  difficult  to  apply  in  its  present  form. 
Also,  several  important  variables  are  not  present.  The  aerodynamic 
forces  depend  not  only  on  the  local  angle  of  attack  of  the  wing  but 
also  on  the  fluid  velocity  and  density.  None  of  these  terms  are 
evident  in  the  equilibrium  expression.  If  the  equilibrium  equation 
is  written: 


cf  -  3  (2) 

where  c  is  the  flexibility  influence  coefficient,  the  displacements  are 
seen  to  rely  upon  the  applied  forces  as  well  as  the  complicated 
structure.  It  is  evident  that  the  equilibrium  equation  must  be 
expanded  and  simplified. 

The  wing  will  be  simplified  by  modeling  it  with  a  finite  element 
model.  ThiB  will  discretize  the  wing  and  allow  the  equilibrium 
equation  to  be  expressed  in  matrix  form.  Then  a  doublet  lattice 
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aerodynamic  method  will  be  applied  to  determine  the  aerodynamic  forces 
as  functions  of  the  local  angles  of  attack  of  the  different  wing 
sections.  This  aerodynamic  model  will  than  be  combined  with  the 
structural  model  so  that  the  eigenvalue  problem  can  be  expressed  by 
a  single  equation.  The  velocity  at  which  the  wing  diverges  will  be 
the  eigenvalue  and  the  rotational  displacements  about  the  y-axis 
(local  change  in  angle  of  attack)  will  become  the  eigenvector. 

Structural  Ming  Model 

The  first  step  in  developing  a  workable  relation  is  to  model 
the  wing  with  a  finite  element  model.  The  wing,  which  has  an  infinite 
number  of  degrees  of  freedom  as  a  continuous  structure,  is  modeled 
as  a  discrete  set  of  points.  At  each  of  these  points,  or  nodes,  the 
structure  is  allowed  to  have  six  degrees  of  freedom,  three  displace¬ 
ments  and  three  rotations.  The  nodes  are  located  at  key  positions  on 
the  wing  so  that  they  can  be  connected  with  simple  structural  elements 
to  approximate  the  actual  structural  members.  Ribs  and  spars  are 
replaced  by  a  combination  of  rod  and  shear  panel  elements.  Skin 
surfaces  are  replaced  by  membrane  or  thin  plate  elements.  In  this 
manner  the  entire  structure  is  broken  down  into  more  workable  elements. 
Figure  1  illustrates  the  finite  element  model  and  also  defines  the 
coordinate  system  used  to  define  the  structure  and  its  degrees  of 
freedom. 

Analytical  methods  are  applied  to  the  individual  elements  to 
distribute  an  effective  stiffness  to  its  nodes.  This  is  done  so  that 
when  the  elements  are  connected,  their  effective  stiffness  terms  add 
to  form  a  stiffness  matrix  for  the  entire  wing  structure.  The 
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stiffness  matrix  for  both  an  element  and  the  entire  structure  can  be 


defined  in  the  following  manner: 

The  elements  k^j  of  the  stiffness  matrix 
are  the  infinitesimal  changes  in  the  forces 
(Fj_)  required  at  the  i1"*1  nodes  to  infinitesimally 
change  the  displacement  at  the  node  j  and  hold 
all  other  displacements  constant,  i.e. 

.  3F* 

The  stiffness  matrix  will  vary  in  shape  and  value  from  structure  to 
structure,  depending  on  the  type  and  number  of  elements  comprising 
the  structure.  However,  the  elemental  stiffness  matrices  for  each 
type  of  element  always  have  the  same  form  and  only  depend  upon  the 
material  and  geometrical  properties  of  the  element.  For  example, 
the  rod  element  has  the  following  stiffness  matrix: 

r.  -i  EA  r-1  n 

L  krod  J  “  L  L  -1  1  J 


The  only  variables  in  thiB  matrix  are  E,  A,  and  L.  They  are  the 
material  and  geometrical  properties  of  the  rod.  Likewise,  shear  panel, 
plate,  and  membrane  element  matrices  can  be  shown  to  depend  upon  only 
geometric  and  material  properties.  (The  elemental  stiffness  matrices 
for  these  elements  will  not  be  reproduced  here  due  to  their  complex 
form.  References  S  and  6  contain  a  more  detailed  account  of  Che 
development  of  the  stiffness  matrices  and  the  finite  element  method 
in  generaL.)  Because  of  this  dependence  on  only  geometric  and  material 
properties,  the  finite  element  method  can  be  readily  programed  into 


a  computer.  There  are  many  widely  used  programs  that  require  only 
the  appropriate  geometric  and  material  properties  of  the  individual 

|  elements  to  assemble  a  stiffness  matrix  for  the  entire  structure. 

I 

Including  the  nodal  forces  might  permit  a  detailed  stress  analysis 
to  be  made.  One  such  program  is  NASTRAN. 

( 

*  NASTRAN  was  developed  by  NASA  to  analyze  virtually  any  structure. 

"  A  wide  variety  of  structural  elements,  ranging  from  the  simple  bar 

element  to  the  complicated  quadrilateral  plate  element,  can  be  used 
to  analyze  both  static  and  dynamic  problems.  It  has  the  capability  to 
find  stresses,  displacements,  natural  modes,  and  natural  frequencies 
and  can  also  be  applied  to  wing  structures  to  calculate  flutter 
speeds.  NASTRAN  is  widely  used  within  the  aerospace  industry  for 
structural  analysis. 

Within  the  framework  of  this  study,  NASTRAN  will  be  used  to 
formulate  the  stiffness  matrix  for  a  test  wing  and  calculate  the 
appropriate  flexibility  matrix.  This  flexibility  matrix  is  the 
inverse  of  the  stiffness  matrix.  It  can  be  found  by  either  inverting 
the  stiffness  matrix  or  applying  the  following  definition: 

Cj  ^  .  .  .  the  i1"*1  column  of  the  flexibility 
matrix  is  a  displacement  configuration  for 

A.  h 

the  j  nodes  resulting  from  a  force  system 
'  of  F^  ■  unity  and  all  other  nodal  forces 

zero . 

Flexibility  influence  coefficients  are  used  quite  often  in  aerodynamic 
analysis  and  will  be  used  to  represent  the  structure  in  this 
development  of  the  eigenvalue  problem. 
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The  finite  element  model  has  effectively  transformed  the 
equilibrium  equation  to  matrix  form: 


<*i>-  [^j]  <V 


(Sj) 

(4) 

<y 

(5) 

or 


The  flexibility  matrix  of  Eq  5  can  be  partitioned  in  such  a  manner 
as  to  separate  the  different  types  of  degrees  of  freedom.  Likewise, 
the  force  and  displacement  vectors  can  be  partitioned  in  a  similar 
manner  to  produce: 
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This  system  can  be  reduced  and  simplified  further  by  making  several 
important  assumptions  regarding  the  structure  and  the  aerodynamic 
forces. 

The  first  assumption  is  that  the  wing  1b  sufficiently  rigid  in 

the  X  -  Y  plane  so  that  the  inplane  displacements  are  small  and  can 

be  neglected.  This  eliminates  the  x,  y,  and  0  degrees  of  freedom. 

z 
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Equation  6  is  reduced  to: 


The  second  assumption  is  that  the  aerodynamic  forces  act  as  a 

distributed  pressure  over  the  upper  surface  of  the  wing.  This  allows 

the  pressure  to  be  represented  by  a  finite  set  of  forces  acting  only 

at  the  upper  nodes ,  eliminating  the  lower  nodes  from  the  equilibrium 

equation.  However,  the  entire  structure  will  still  be  represented  by 

the  values  of  the  remaining  flexibility  coefficients.  The  remaining 

forces  in  the  x  and  y  directions  have  already  been  eliminated  with 

the  previous  assumption  leaving  only  a  set  of  aerodynamic  forces  In 

the  z  direction  acting  on  the  nodes  of  the  wing's  upper  surface.  Thane 

forces  may  induce  a  resultant  moment  about  either  the  X  or  Y  axis  of 

the  wing,  however;  there  are  no  aerodynamic  moments  applied  at  the 

individual  nodes  (i.e.  M  ■  M  ■  0).  This  results  in  the  following 

X  y 

set  of  3N  (N  is  the  number  of  upper  surface  nodes)  simultaneous 
equations: 


3NxN  Nxl  3Nxl 
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Note,  however,  that  this  set  of  equations  can  be  separated  Into  three 
Independent  sets  of  N  simultaneous  equations, 


£  Czz  ^NxN  i 

\ 

tFzW 

■ 
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* 0y*Nxi 

(9c) 

\ 

Each  of  these  can  bdj,  solved  independently  from  the  others. 

'  l 
\\ 

Equation  9c  is  ,bf  most  importance  to  the  divergence  problem.  As 
will  be  shown  next,  the  aerodynamic  forces  are  related  to  the  local 
angle  of  attack  of  the  wing.  This  local  angle  of  attack  can  be 

related  to  the  displacement  6y  either  directly  or  indirectly, 

I 

depending  upon  the  njbnner  in  which  the  aerodynamic  and  structural 

> 

models  are  combined,!  Consequently,  Eq  9c  is  most  suitable  for 
developing  into  the  eigenvalue  problem. 


Aerodynamic  Model 


The  next  step  In  developing  the  eigenvalue  problem  la  to  develop 
a  workable  aerodynamic  model.  The  foundation  for  this  model  will  be 
the  tangential  flow  condition  aa  developed  by  Aahley  and  Landahl 
(Ref  1) .  It  relates  the  local  angle  of  attapk  to  the  vortex 
distribution  which  in  turn  can  be  related  to  the  lift  distribution. 
However,  this  relation  is  difficult  to  work  ^ith  and  will  be  simplified. 
The  simplification  ia  accomplished  with  a  doublet  lattice  modeling 
technique.  This  will  leave  the  tangential  flow  condition  in  a  form 
that  can  be  combined  with  Eq  9c  to  form  the  'eigenvalue  problem.  The 

V 

t 

doublet  lattice  method  ia  described  in  more  detail  in  Appendix  A. 

i 

The  wing  is  regarded  as  a  thin  wing.  This  allows  it  to  be 

> 

modeled  aa  a  sheet  of  pressure  doublets.  Ttie  strength  of  a  pressure 

"l 

doublet  corresponds  to  the  strength  of  a  hotjseahoa  vortex  as  It  is 
shrunk  to  a  point.  This  includes  the  effects  of  the  trailing  portions 
of  the  vortex  and  allows  the  tangential  flo^  relation  to  be  formulated 
(Ref  I) : 

1  f  ( 

v(K,y,0)-j-j  j 
^wing 

where  Y(xi»y],)  the  strength  of  the  pressure  doubi.'t  or  vortex 
sheet  at  the  point  (x^,y^).  This  equation  Relates  the  component  of 

i 

the  fluid  velocity  in  the  z  direction  at  this  wing's  surface  to  the 
pressure  doublet  distribution.  For  wings  at  small  angles  of  attack, 
the  w  component  of  the  fluid  velocity  can  be  approximated  by: 


i  + 


(x-xj,)2  d  (y-y^)' 


dxrdy. 


f'  n\ 
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w(x,y,0)  ■  U.  -U„a(x,y)  (11) 

This  allow*  tha  local  angle  of  attack  of  a  point  (x^y^O)  to  b* 
written  in  terms  of  tha  doublet  distribution: 


<*(x1,yi> 


r  f  y(*i»yi> 
4irU.  J  J  (y^yj.)2 

swing 


(Xi-Xj) 

/Oti-xjT2  +  (yi-yi): 


dx^yx  (12) 


For  the  divergence  problem,  it  is  advantageous  to  obtain  tha 
vox  tax  strength  (r)  as  a  function  of  the  wing's  geometry  and  flow 
conditions,  The  reason  for  this  is  that  the  vortex  distribution  is 
related  to  the  normal  lift  distribution  by:  (Ref  A) 

L  ■  J  P»  u«.  Y(*.y>  dx  (13) 

o 

where  L  is  the  lift  per  unit  span.  This  relates  the  lift  to  the 
wing's  geometry.  If  discretized  properly,  the  lift  forces  can 
be  combined  with  the  equilibrium  equation  to  form  the  eigenvalue 
problem.  However,  the  above  tangential  flow  integral  equation  is 
very  difficult  to  invert  using  analytical  methods.  If  the  wing  is 
modoled  with  the  doublet  lattices  method,  the  integral  equation  is 
transformed  to  a  matrix  equation.  The  vortex  or  doublet  distribution 
can  then  be  found  by  simply  inverting  a  matrix  of  aerodynamic 
coef f icientb. 

The  wing  is  divided  into  a  finite  number  of  aerodynamic  panels 
(see  ^ig  2).  Each  panel  cm  'ins  a  control  point  (x^,y^),  at  which 


the  angle  of  attack  is  known,  and  a  point  (Xj.yj)  for  locating  a  vortex 
singularity.  For  an  individual  panel  the  vorticity  (circulation  per 
unit  area)  is  approximated  by  a  circulation  per  unit  span  acting  along 
a  constant  cord  line  through  the  point  (Xj,yj).  This  is  achieved  by 
representing  Y(x^,y^)  with  a  Dirac  delta  function 

Y^xl,yl^  "  r(yi*  6<xi  “  xj*  f1*) 

Substituting  this  into  equation  12  and  noting  that 

B 

f(t)  $<t  -  tQ)  dt  -  f(t0)  (15) 

A 

for  A  <  tQ  <  B 

allows  the  integration  over  x^  to  be  performed: 


u«.  “  7 


i  N  ( 

57  J-l  ^ 


yrsj 


r(yl) 

(yi‘yi)2 


l  + 


(xi'xj)2+(yi-yi)2J 


dyi 


N  -  number  of  panels 


(16) 


The  circulation  r(y^)  can  be  brought  outside  the  integral  if  it  is 
assumed  to  be  «  nnm  .uii  across  the  span  of  each  panel.  The  remaining 
integral  is  an  improper  integral  oi  the  second  kind.  For  a  number  of 
the  panela  there  are  singularity  points  at  y^  *  y^  However,  it  can 
be  integrated  using  Mangler's  principal-value  technique  (Ref  1). 


15 


{ r . }  *  the  strengths  of  circulation  singularities 
^  located  at  the  points  (Xj,yj)  j»l,  .  .  .  N 

The  circulation  distribution  is  obtained  by  simply  multiplying 
the  equation  by  the  inverse  of  the  aerodynamic  matrix: 

<y  “  u»  £>i,j  r1  (y  as) 

For  an  individual  panel,  the  lift  per  unit  span  is  related  to  the 
circulation  by: 

Lj  -  P«  rj  (20)' 

or  in  matrix  form: 

<>•! )  -  p. ».  o-j )  ■  p.  ".2  [  Vj  r1  v 

The  total  lift  of  each  panel  is  found  by  multiplying  Lj  by  the  span 
length  of  the  panel. 

<V  -  P„  U,»  [  2sj  ‘J  -  P.  U.2  I.28J  ]  [  A1(j  Yl  {at}  (22a) 

If  the  span  length  1b  the  same  for  all  panels: 

{y  ■  2pw  S  <y  -  2P»  "J  s  [  'J'1  {at}  (22b) 

This  effectively  replaces  the  continuous  pressure  distribution  with 
a  set  of  forces  acting  at  the  singularity  point  locations.  Each  force 
depends  upon  the  angle  of  attack  at  the  control  points  of  the  other 
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panels,  which  may  vary  from  panel  to  panel.  After  the  angles  of 
attack  of  the  aerodynamic  panels  are  related  to  the  rotational 
degrees  of  freedom  at  the  structural  nodes,  the  aerodynamic  and 
structural  models  can  be  combined  to  form  the  eigenvalue  problem. 

Combination  of  Aerodynamic  and  Structural  Models 

The  final  obstacle  in  formulating  the  eigenvalue  problem  is 
encountered  when  combining  the  aerodynamic  and  structural  models. 

The  local  angle  of  attack  at  the  aerodynamic  panel  control  points 
must  be  related  to  the  rotational  degrees  of  freedom  at  the 
structural  nodes.  It  may  be  helpful  to  consider  a  rectangular  wing 
which  has  a  simple  flat  plate  for  its  load  bearing  structure.  The 
structure  is  modeled  by  a  set  of  nodes  connected  by  quadrilateral 
plate  elements  (see  Fig  3).  An  aerodynamic  planform  can  be 
overlaid  upon  the  structural  model  such  that  the  aerodynamic 
singularity  points  coincide  with  the  structural  nodes.  This  allows 
the  aerodynamic  forces  to  be  applied  directly  to  the  nodes. 

There  are  several  ways  in  which  the  local  angle  of  attack  of  the 
control  points  can  be  related  to  the  rotational  degrees  of  freedom. 

One  method  Is  to  assume  that  each  aerodynamic  panel  acts  as  a  flat 
plate.  This  illows  the  angle  of  attack  at  a  control  point  to  be 
approximated  by  the  angle  of  attack  at  the  singularity  point  of  the 
same  panel.  Since  the  singularity  point  coincides  with  the  node,  the 
angle  of  attack  at  the  singularity  point  is  the  rotation  at  the  node. 
Thus,  the  aerodynamic  and  structural  models  could  be  combined  directly 
to  form  the  eigenvalue  problem.  However,  this  method  may  incur 
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problems  because  it  allows  discontinuities  across  the  aerodynamic 
panel  borders. 

A  second  method  assumes  that  the  aerodynamic  planfora  is  rigid 
in  the  X  direction.  The  angle  of  attack  for  a  given  span  station  is 
found  by  dividing  the  difference  between  the  displacements  of  the 
leading  edge  and  trilling  edge  structural  nodes  by  the  structural 
chord  length.  This  angle  is  used  for  the  angle  of  attack  of  all 
control  points  along  the  chord  of  that  particular  span-wise  station. 

In  this  manner  the  aerodynamic  model  is  indirectly  combined  with  the 
structural  model  to  form  the  eigenvalue  problem. 

Both  of  the  above  methods  provide  similar  resultB.  They  use  an 
effective  angle  to  approximate  both  the  local  svgl...  of  attack  at  the 
aerodynamic  control  points  and  the  rotational  degrees  of  freedom  at 
the  structural  nodes.  The  final  form  of  the  eigenvalue  problem  is 
formulated  by  substituting  the  lift  forceB  of  Eq  22b  into  the 
equilibrium  Eq  9c,  and  replacing  both  the  local  angle  of  attack  a 
and  the  rotational  degrees  of  freedom  6y  with  the  effective  angle 

9eff 

A  (S)  (-^-)  [C0z]  [A^  J-1  <Geffl)  -  {9effj}  (23) 


The  lowest  value  of  the  eigenvalue  (1/2  is  the  divergence 

dynamic  pressure,  from  which  the  divergence  speed  can  be  readily 
determined.  The  flexibility  and  aerodynamic  influence  coefficients 
are  easily  determined  with  the  help  of  a  computer,  making  the 
eigenvalue  problem  particular  suitable  to  computer  aided  eigenvalue 
analysis  techniques. 
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t  A  wing  with  an  aluminum  box  structure  (complex  structures  are 

}  ‘ 

discussed  in  the  next  section)  was  analyzed  using  both  methods  for 
'  obtaining  the  effective  angle  6eff  Figure  4  compares  these  results 

with  the  results  obtained  by  using  NASTRAN1 a  capabilities  to  determine 
the  divergence  speed.  (The  NASTRAN  computer  program  is  discussed  in 
more  detail  later.)  As  is  evident  from  Fig  4,  both  methods  for 
combining  the  aerodynamic  and  structural  models  produce  nearly 
identical  results.  The  variation  of  divergence  speed  with  sweep  angle 
also  correlates  with  the  results  obtained  from  NASTRAN  if  the  differ¬ 
ences  in  the  two  methods  and  NASTRAN  are  accounted  for  (NASTRAN  used  a 
strip  theory  model  and  assumed  a  symmetric  lift  distribution,  resulting 
in  lower  values),  indicating  the  method  to  be  valid.  The  second  method 
for  approximating  an  effective  angle  of  attack  will  be  used  in  the  anal- 
i"  yBis  of  the  test  wing.  It  will  provide  an  accurate  approximation  for 

the  local  angle  of  attack  since  the  test  wing  has  a  simple  box  structure. 

Effects  of  Complex  Wing  Structures 

For  a  simple  rectangular  wing,  Buch  as  the  flat  plate  example  (Fig 
3),  it  is  easy  to  visualize  how  the  aerodynamic  grid  might  be  directly 
overlaid  on  to  the  structural  nodes.  However,  sweeping  the  wing  forward 
complicates  the  problem  since  it  can  be  done  in  several  ways. 

One  method  is  to  "sweep"  the  wing  forward  in  such  a  manner  that 
the  root  and  the  tip  remain  parallel  to  the  flow  (Fig  5).  This 
method  allows  the  aerodynamic  representation  to  remain  simple  but 
"skews"  the  structural  elements  so  that  the  nodes  remain  aligned  with 
the  flow.  A  second  method  is  lo  "rotate"  the  wing  forward  (Fig  6). 

This  allows  the  structural  representation  to  stay  the  same  but  "skews" 
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—  Method  1  (Rotations  of  Nodes) 

. . . -  Method  2  (Difference  of 

Displacements  of 
Nodes) 

-  NASTRAN  (Strip  Theory) 


— i - 1  -  I  1  ■■■"■"I  1 . "i -  '1  .  |  ■  f 

10  20  30  40 

Forward  Sweep  Angle  (-A)0 

Figure  4.  Comparison  of  Methods  for  Obtaining 
Divergence  Speeds. 
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the  aerodynamic  boxes  so  that  the  aerodynamic  grid  points  coincide 
with  the  structural  nodes.  With  this  method,  the  aerodynamic  panels 
do  not  align  in  the  chord-wise  direction.  This  may  or  may  not  be 
troublesome,  depending  on  the  manner  in  which  the  effective  angle  of 
attack  is  found. 

The  flat  plate  wing  of  Fig  3  is  not  affected  by  the  pitfalls 
of  either  method.  It  is  represented  by  in-plane  plate  elements  which 
do  not  have  a  strong  dependence  upon  their  shape  to  provide  out-of¬ 
plane  bending  stiffness.  Unfortunately,  this  simple  wing  model  is 
not  very  prevalent  among  airplanes.  A  more  common  wing  structure 
might  be  represented  by  the  box  structure  shown  in  Fig  1.  This 
structure  is  affected  by  the  manner  in  which  it  is  swept  forward. 
Requiring  the  nodes  to  be  aligned  with  the  flow  as  the  wing  is  swept 
forward,  also  aligns  the  chord -wise  shear  panels  with  the  flow.  This 
changes  the  stiffness  characteristics  of  the  structure  as  it  is  swept 
forward.  As  compared  to  the  rotated  wing,  the  swept  wing  will  have 
a  larger  torsional  stiffness  and  a  smaller  bending  stiffness.  However, 
for  the  simple  box  wing  to  be  studied,  this  difference  will  be  small. 
Thus,  the  first  method,  in  which  the  wing  is  "swept"  forward,  will  be 
used  to  analyze  the  test  wing. 

Remarks 

The  method  described  above  for  determining  the  divergence  speed 
may  appear  to  be  contrived.  The  method  might  not  be  applicable  to 
wings  of  arbitrary  shape  and  the  test  wing  might  seem  to  be  designed 
around  this  method.  However,  it  should  be  pointed  out  that  the  purpose 
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of  this  study  is  to  quantitatively  determine  how  the  divergence  speed 
of  a  wing  at  a  given  sweep  angle  is  affected  by  the  orientation  of  the 
composite  laminates  making  up  the  wing's  structure.  The  material  and 
structural  properties  of  the  forward  swept,  composite  wing  are 
represented  by  the  flexibility  matrix.  The  aerodynamic  properties 
are  represented  by  the  aerodynamic  matrix.  Any  variation  in  the 
structure  or  sweep  angle  is  reflected  in  these  two  matrices,  causing 
the  divergence  speed  represented  by  the  eigenvalue  to  vary  accordingly. 

This  method  can  be  applied  to  wings  in  either  symmetric  or 
antisymmetric  flight.  Symmetry  is  accounted  for  by  adding  the  effects 
of  the  left  wing  to  the  aerodynamic  matrix.  The  divergence  speed  for 
a  wing  in  symmetric  flight  will  then  be  less  than  that  for  the  same 
wing  with  an  antisymmetric  lift  distribution.  However,  this  has  no 
affect  on  the  manner  in  which  composite  fiber  orientation  affects  the 
divergence  speed.  Therefore,  an  antisymmetric  lift  distribution  is 
assumed  because  it  simplifies  the  formulation  of  the  aerodynamic  matrix. 
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III.  Analysis  of  Example  Wings 


Introduction 

The  previous  chapter  developed  an  eigenvalue  problem  for  finding 
the  divergence  spaed  of  forward  swept  wings.  A  computer  program  for 
formulating  and  solving  this  eigenvalue  problem  was  written  and 
applied  to  a  series  of  example  wings.  The  computer  program  used 
NASTRAN  to  determine  the  flexibility  matrix;  calculated  the 
aerodynamic  matrix  from  Eq  18;  solved  the  eigenvalue  problem  with 
an  inverse  power  solution  technique;  and  calculated  the  divergence 
speed  from  the  lowest  eigenvalue. 

First,  an  aluminum  wing  was  analyzed.  Its  divergence  speed  was 
calculated  for  different  sweep  angles  to  provide  a  basis  of  comparison 
with  other  methods.  Than,  the  upper  and  lower  surfaces  of  the  wing 
were  changed  to  composite  laminates.  For  a  given  sweep  angle  the 
composite  fiber  orlentationu  were  varied.  The  divergence  speeds 
were  calculated  for  the  various  fiber  orientations  and  compared.  The 
results  of  this  analysis  are  presented  in  Chapter  5.  A  more  detailed 
description  of  this  procedure  follows. 

Description  of  Example  Wings 

The  wings  to  be  analyzed  have  basically  a  cantilevered  box 
structure,  as  is  shown  in  Fig  7.  They  are  composed  of  seven  ribs 
and  two  spars  each  and  are  encased  within  a  thin  airfoil,  The 
particular  dimensions  of  both  the  aerodynamic  planform  and  the 
structural  box  are  shown  in  the  figure.  The  wings  are  swept  forward 
such  that  the  total  area  of  the  wing  remains  constant.  The  lengths 
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of  the  spars  remain  constant  while  the  rib  lengths  are  allowed  to 
increase  as  the  wings  are  swept  forward.  The  ribB  alwaya  remain 
parallel  to  the  flow.  This  allows  the  basic  structural  box 
dimensions  to  remain  constant  as  the  wings  are  swept  forward.  It 
also  permits  the  structure  and  aerodynamic  planform  of  each  wing  to 
be  modeled  in  a  manner  that  provides  for  the  easy  combination  of  the 
two. 

Structural  Model 

The  wings  defined  above  are  each  represented  by  28  nodes,  14  on 
the  upper  surface  and  14  directly  below  them  on  the  lower  surface. 

Each  node  is  located  at  the  junction  of  a  rib,  a  spar,  and  the  upper 
or  lower  surface.  The  four  nodes  at  the  root  of  each  wing  are  totally 
fixed  to  cantilever  the  structure.  The  rib  and  spar  members  are 
represented  by  shear  panels  and  rod  elements.  The  rod  elements  are 
distributed  so  that  each  shear  panel  has  rod  elements  on  all  four 
aides  (some  of  the  rod  elements  are  shared  between  shear  panels) . 

These  el». monts  have  aluminum  material  properties.  The  material 
properties  and  the  cross  sectional  dimensions  of  these  elements 
remain  -he  same  for  every  case  considered. 

The  upper  and  lower  surfaces  are  represented  by  thin  plate 
elements •  For  the  aluminum  wings,  there  are  12  plate  elements, 

6  ou  the  upper  surface  and  6  on  the  lower  surface.  Each  element  Is 
0325  in.  thick.  For  the  composite  wings,  there  are  120  plate 
elements.  Each  aluminum  plate  element  is  replaced  by  10  composite 
plate  elements  .00525  in.  thick.  Each  element  represents  a  composite 
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l.mina  and  the  10  elements  are  connected  to  the  same  4  nodes.  This 
builds  up  a  10  layer  composite  laminate  for  representing  the  wing 

Si 

surface.  The  principal  fiber  directions  of  each  of  these  layers  can 
be  oriented  differently  to  vary  the  structural  characteristics  of  the 
wing. 

Two  cases  will  be  considered.  These  cases  ere  defined  in  Fig  8 
and  Fig  9.  For  both  cases ,  the  principal  fiber  direction  6  will  be 
varied  from  0°  to  50°  ahead  of  the  structural  sweep  axis  for  several 
wing  sweep  angles.  The  computer  program  NASTRAN  is  used  to  calculate 
the  flexibility  matrix  coefficients  for  each  instance.  The  geometric 
and  material  properties  of  the  structure  and  its  elements  ere  input  to 
NASTRAN  via  the  appropriate  data  carda  to  permit  NASTRAN  to  formulate 
the  elemental  stiffness  matrices  and  assemble  them  into  the  structural 
1  stiffness  matrix. 

The  different  methods  for  combining  the  structural  and  aerodynamic 
models  require  different  flexibility  matrices.  For  both  methods,  a 
unit  load  is  applied  to  the  1^  upper  surface  node,  while  all  other 
forces  are  zero.  This  is  done  in  turn  to  each  of  the  upper  surface 
nodes.  For  each  load  case  NASTRAN  solves  the  equilibrium  equations: 

{F}  ■  [  k  ]  {6}  (4) 

for  all  nodal  displacements.  Only  the  displacements  for  the  upper 
surface  nodes  are  output. 

The  direct  combination  method  uses  the  rotations  output  by  NASTRAN 
.  for  the  flexibility  matrix.  The  rotations  about  the  y-axis  of  the 

upper  surface  nodes,  caused  by  a  unit  load  at  the  i^  node, 


29 


Nodes 

(28) 


Rod  Elements  - 
(on  edges  of 
shear  panels) 
(48) 


'Composite  Fibers 
1002  along  0 


Thin  Plate 
Elements 
(12  aluminum 
or  120 
composite) 


Composite  Fibers 
702  along  0 
102  along  45 
102  along  90 
102  along  135 


Figure  9.  Case  2;  702  of  Fibers  Oriented  Along  the  0  Direction 
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corresponds  to  the  element  of  the  flexibility  matrix.  Thus, 

NASTRAN  directly  computes  an  entire  column  of  the  flexibility  matrix 
for  each  unit  load  force  set. 

The  indirect  combination  method  requires  additional  calculations. 
For  a  unit  load  at  the  it*'  node,  the  difference  in  z  displacements 
between  the  J  and  j  +  1  nodes  (The  j  and  j  +  1  nodes  are,  respectively, 
the  leading  and  trailing  edge  nodes  of  a  given  span-wise  station.)  is 
divided  by  the  structural  distance  between  the  j  and  j  -I-  1  nodes. 

This  calculated  angle  is  then  used  for  the  and  i  elements 

of  the  flexibility  matrix. 

The  Indirect  method  is  used  for  the  eigenvalue  analysis  of  the 
various  wings. 


Aerodynamic  Model 

The  aerodynamic  planform  is  attached  to  the  upper  surface  of  the 
structure  and  is  sectioned  into  symmetric  panels  as  shown  in  Fig  10. 

The  vortex  singularities  are  located  at  the  quarter-chord  mid- 
span  point  of  each  panel.  This  allows  them  to  coincide  with  the 
upper  surface  nodes.  The  control  points  are  located  at  the  three- 
quarter  chord  mid-span  position  of  each  panel.  The  computer  program 
automatically  calculates  these  points  for  a  given  sweep  angle.  It 
then  uses  these  coordinates  to  calculate  the  aerodynamic  matrix  elements 
according  to  Eq  17  and  18.  The  flexibility  matrix  calculated  by 
NASTRAN  is  input  into  the  program  and  multiplied  by  the  inverse  of 
the  aerodynamic  matrix.  This  sets  up  the  eigenvalue  problem: 


4S  (  PlpU°°  )  f  CA  ]  (0)  -  {0} 
2 
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Solution  of  the  Eigenvalue  Problem 

An  inverse  power  method  is  used  to  solve  the  eigenvalue  problem. 
The  method  described  by  Hornbeck  (Ref  3)  was  coded  into  FORTRAN  and 
applied  to  the  product  matrix  calculated  above.  This  method  is  an 
Iterative  process  that  converges  to  the  lowest  eigenvalue. 

The  first  step  is  to  guess  a  solution  for  the  eigenvector.  This 
guess  is  thfn  multiplied  by  the  product  matrix.  The  resultant  vector 
is  normalized  to  the  value  of  the  last  element  of  its  vector.  This 
normalized  vector  is  then  multiplied  by  the  produc  matrix.  The 
resultant  vector  is  then  normalized  and  again  multiplied  times  the 
product  matrix.  This  process  continues  until  the  normalized  vectors 
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converge  Co  the  eigenvector.  Similarly ,  Che  normalizing  factor  will 
converge  to  the  inverse  of  the  eigenvector  (  y  ) .  This  method  always 
converges  to  the  lowest  eigenvalue,  if  it  exists,  usually  in  less  than 
6  iterations.  From  Eq  24  the  eigenvalue  is: 

PoaUoo^ 

X  -  4S  (  — —  )  (25) 

The  divergence  speed  is  found  by  solving  Eq  25  for 

UD1V  “  U«  “  r.7sp^"|1/2  <26> 

A  listing  of  both  the  NASTRAN  data  deck  and  the  program  written  to 
calculate  the  aerodynamic  matrix  and  solve  the  eigenvalue  problem  is 
included  in  appendix  B.  The  results  of  this  analysis  are  presented 
in  the  next  chapter. 
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IV.  Results 


The  wing  defined  in  the  previous  chapter  is  used  to  study  how 
the  divergence  speed  varies  with  composite  fiber  orientation.  First, 
divergence  speeds  of  the  aluminum  wing  at  sweep  angi.es  of  0°,  15°, 

25°,  35°,  and  45°  are  calculated  using  the  computer  program  described 
above.  These  values  are  plotted  in  Fig  4  and  are  used  as  a  basis  for 
comparison  with  other  methods.  For  the  two  composite  wing  cases 
(Fig  8  and  9),  the  fiber  orientation  (6)  is  varied  in  5°  increments 
from  0°  to  50°  ahead  of  the  structural  sweep  angle.  The  divergence 
speed  for  each  increment  is  computed  and  compared  to  the  divergence 
speed  for  the  composite  wing  with  the  fibers  oriented  at  6  ■  0.  This 

nondimenslonalizes  the  divergence  speed  with  the  ratio  ,  6 —  . 

UDIV  0-0 

Figures  11  and  12  plot  this  divergence  speed  ratio  vereuB  the  composite 
orientation  angle  for  the  five  sweep  angles.  Figure  11  considers  the 
case  where  100%  of  the  composite  fibers  are  oriented  in  the  0  direction 
and  Fig  12  considers  the  case  where  only  70%  of  the  fibers  are 
oriented  in  the  8  direction  (the  remaining  fibers  are  oriented  at  90° 
and  +  45°  to  the  structural  sweep  angle) .  Figure  13  shows  how  UDIV 
varies  with  sweep  angle  for  6  -  0. 

From  these  results  it  is  evident  that  the  divergence  problem  of 
forward  swept  wings  can  be  effectively  controlled  with  laminated 
composites.  As  the  principal  fiber  direction  is  swept  ahead  of  the 
wing  sweep  angle,  the  divergence  speed  is  seen  to  lncreeae  and  reach 
a  maximum.  For  wings  of  lesser  sweep  angles,  the  divergence  problem 
is  eliminated  completely.  This  is  true  for  both  cases,  although, 


34 


in  varying  degrees. 

With  the  optimal  choice  of  fiber  orientation  angle,  the 
divergence  speeds  for  case  1  (Fig  11)  are  increased  by  a  factor 
of  more  than  1.5  for  a  45°  swept  wing,  and  totally  eliminated  for 
wings  swept  25°  or  less.  Optimal  fiber  orientations  are  between 
15°  and  20°  for  wings  swept  more  than  15°.  The  divergence  speeds 
for  case  2  (Fig  12)  are  not  affected  as  drastically  as  they  wero 
for  case  1,  however;  the  same  trends  were  obtained.  Maximum  values 
for  case  2  for  wings  swept  15°  or  more  occur  for  principal  fiber 
orientations  between  15°  and  25°. 

In  both  cases,  the  divergence  speed  varied  unexpectedly  for 
wings  of  small  sweep  angles.  Unswept  wings  encountered  two  different 
maximums  in  the  divergence  speed  variation  at  two  different  fibar 
orientations,  around  G  -  10°  and  G  ■  40°.  No  attempt  has  been  made 
to  explain  this  phenomenon  other  than  to  suggest  that  it  might  be 
caused  by  the  coupling  of  the  torsion  and  bending  modes  of  divergence. 


V.  Conclusion 


Summary 

An  algenvalua  problem  was  formulated  to  study  the  effects  of 
composite  material  fiber  orientations  on  the  divergence  speed  of 
forward  swept  wings.  This  was  accomplished  by  relating  the  aerodynamic 
forces  to  the  wing  deformations  with  the  tangential  flow  condition. 
These  deformations  were  related  to  structural  properties  by  applying 
the  static  equilibrium  equation.  The  aerodynamic  and  structural 
properties  were  modeled  with  a  doublet  lattice  method  and  a  finite 
element  method,  respectively.  These  modeling  techniques  simplify  the 
appropriate  force-displacement  relations  and  allow  the  relations  to  be 
written  in  matrix  form.  The  two  force-displacement  relations  ware 
combined  to  form  the  eigenvalue  problem  of  Eq  23. 

A  computer  program  was  written  which  used  NASTRAN  to  calculate  a 
flexibility  matrix  for  the  structure.  It  calculated  the  aerodynamic 
matrix  from  the  doublet  lattice  method  and  solved  the  eigenvalue 
problem  with  an  inverse  power  technique.  The  program  was  applied 
to  two  cases  of  swept  wings.  For  the  first  case,  100%  of  the  composite 
fibers  were  varied  from  0°  to  50°  ahead  of  the  structural  sweep  axis. 
For  the  second  case,  only  70%  of  the  compooite  fibers  ware  varied 
while  the  remaining  fibers  were  held  at  90°  and  +  45°  to  the 
structure-swept  axis.  Sweep  angles  from  0°  to  45°  were  studied. 

L'he  calculated  divergence  speeds  for  both  cases  are  compared  in 
Fig  11  and  Fig  12. 
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Conclusions 

Figures  11  and  12  indicate  that  the  divergence  problem  of  foi  -r<: 
swept  winga  can  be  effectively  controlled  and  possibly  eliminated  by 
using  the  proper  composite  fiber  orientations.  The  divergence  speed 
was  found  to  Increase  and  reach  a  maximum  as  the  principal  fiber 
direction  was  swept  forward.  The  maximum  values  occurred  for  fibers 
oriented  between  15°  and  20°  ahead  of  the  structural  sweep  axis  for 
wings  swept  forward  25°  or  more.  This  was  found  to  be  true  for  both 
cases,  although,  to  a  lesser  extent  for  the  case  where  only  70% 
of  the  fibers  were  varied. 

The  results  obtained  by  this  study  concur  with  the  results  of 
both  Krone  (Ref  6)  and  Weisshaar  (Ref  7)  in  affirming  the  feasibility 
of  using  laminated  composites  for  controlling  the  divergence  speed  of 
forward  swept  wings.  However,  no  attempt  has  been  made  to  study  the 
affects  of  composite  materials  on  flutter,  strength,  fatigue  or  other 
problems  of  wing  structures.  Further  study  is  required  in  these 


areas. 
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Appendix  A 

Doublet  Lattice  Aerodynamic  Model 


The  doublet  lattice  method  Incorporates  the  usual  assumptions 
regarding  thin  airfoils  and  lrrotatlonality  to  form  a  tangential  flow 
equation.  This  equation  relates  the  angle  of  attack  of  a  given  point 
on  the  wing  surface  to  the  strengths  of  the  circulation  vortices 
representing  the  wing.  In  most  instances,  the  geometry  of  the  wing  Is 
known  and  the  problem  is  to  determine  the  circulation  distribution. 
This  is  most  efficiently  accomplished  using  matrix  methods.  The  wing 
is  dlscretised  to  allow  the  tangential  flow  equation  to  be  written  in 
matrix  form.  The  equation  can  then  be  inverted  and  the  circulation  at 
a  given  point  found  in  terms  of  the  wing's  geometry  and  the  fluid 
velocity.  From  elementary  aerodynamics,  the  circulation  is  related 
to  the  normal  lift  forces  providing  a  useful  method  of  detarminlng  the 
lift  and  drag  forces  and  the  resultant  moments  acting  upon  a  wing. 

The  wing  is  regarded  as  a  thin  wing.  This  allows  it  to  be 
modeled  as  a  sheet  of  pressure  doublets.  The  strength  of  a  pressure 
doublet  corresponds  to  the  strength  of  a  horseshoe  vortex  as  it  is 
shrunk  to  a  point.  This  includes  the  effects  of  the  trailing  portions 
of  the  vortex  and  allows  the  tangential  flow  relation  to  be  formulated 
(Ref  1): 


,  i  Cf  vUi.yp 

»C«.y.o)-5vJj^yr- 


1  + 


(x-Xj^ 


f. 


(x-Xj^)  -+-<y-y1>  ‘ 


dx1d>,i 


<A-1) 


where  7(x^,y^)  is  the  strength  of  the  pressure  doublet  or  vortex 
aheet  at  the  point  (x^,y^).  This  equation  relates  the  component  of 
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th*  fluid  velocity  in  the  z  direction  et  the  wing's  surface  to  the 
pressure  doublet  distribution.  For  wings  at  small  angles  of  attack, 
the  w  component  of  the  fluid  velocity  can  be  approximated  byt 

5* 

w(x,y,0)  *  U  — *  ■ -U„a(x,y)  <A-2) 

*  Ox 


This  allows  the  local  angle  of  attack  of  a  point  (x^,y^,0)  to  be 
written  in  terms  of  the  circulation  distribution: 


However,  in  most  instances,  the  vortex  distribution  is  not  known 
end  it  would  be  desirable  to  obtain  the  vortex  strength  at  the 
point  (x^y^)  in  terms  of  the  wing  geometry.  This  is  most  effectively 
accomplished  using  matrix  methods.  Therefore,  the  next  step  le  to 
dlscretlse  the  wing  to  permit  the  tangential  flow  condition  to  be 
written  in  matrix  form. 

The  wing  is  divided  into  a  finite  number  of  aerodynamic  panels 
(sea  Fig  2).  Each  panel  contains  a  control  point  (x^y^),  at  which 
the  angle  of  attack  is  known,  end  a  point  <Xj,yj)  for  locating  a 
vortex  singularity.  For  an  individual  panel  the  vortlcity  (circulation 
per  unit  area)  is  approximated  by  a  circulation  per  unit  span  acting 
along  a  constant  chord  lina  through  the  point  < x j , y j > .  This  is 
echisvsd  by  rspresenting  y(x^,y^)  with  s  Dirsc  delte  function: 

m  rCy^  5(xi  -  xj)  (a-4) 
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1 


Substituting  this  into  Eq  A- 3  and  noting  that: 


f(t>  i(t-to)  dt  -  f(t0) 


(A-5) 


for  A  <  tQ  <  B 


allows  tha  Integration  over  x^  to  be  performed: 


yJ+SJ 


U  a 


.1  M  f  J  J  r(y.) 

tx^i)  -  ^  l  f  - - 

yryl 


1  + 


(x1~x1) 


(*1-«J)2+(y1-y1)2 


dVi 


N  *  number  of  panels  (A-6) 


The  circulation  r(y^)  can  be  brought  outside  the  integral  if  it  is 
assumed  to  be  constant  across  the  span  of  each  panol.  The  remaining 
integral  is  an  improper  integral  of  the  second  kind.  For  a  number  of 
the  boxes,  there  are  singularity  points  at  y^  ■  y;, .  However,  it  can 
be  integrated  using  Mangier' s  principal-value  technique.  (Ref  1) 
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u» 


N  rYsj 


v  -s  (yrY' 

yj  j 


1  +  ^rY _ 

,A(x1-xJ)2+(y1-y1) 


j-1  Airfx^Xj) 


Uj-x.)  +  /(xt-x,)2  +  (y^-y.-S .  )2 


yi  -  yi  ‘  sj 


(xi-x.)  +  "(x^x,)2  +  (yi-y,+S  )2 


yt  “  yj  +  Sj 


(A-7) 


This  equation  can  now  be  written  in  matrix  form; 


u~  'V  •  CAi,jl  <V 


U-a> 


where 


-  free  stream  velocity 

{a^}  ■  the  angles  of  attack  of  the  set  of  control 
points  (x^,7^),  i-1,  .  .  .  N 

[At  j]  ■  matrix  of  aerodynamic  influence  coefficients 

A  .  li_  fx  +  yx2  +  (y  -  SQ2 
1,3  4itx  (,  y  -  S. 


x  -t-  x2  +  (y  +  S.)2 


y  +  S, 


y  "  yi  "  yJ 
x  -  ~  *j 


45 


(T. }  ■  the  strengths  of  circulation  oingularities 
J  located  at  the  points  (x^.y^)  J-l,  «  .  .  N 

The  circulation  distribution  is  obtained  by  simply  multiplying 
Eq  A-fl  by  the  inverse  of  the  aerodynamic  matrix: 

{V  "  u«  £>i,j  yl  {oi}  (A“9) 

For  an  individual  panel,  the  lift  per  unit  span  is  related  to  the 
circulation  by; 

-  p. ». <*-1« 

or  in  matrix  form: 

V  ■ ». »y  -  <>.  ».2  [  j'1  <v  <*-u> 

This  provides  the  lift  per  unit  span  for  the  individual  aerodynamic 
panels.  The  total  lift  on  each  panel  is  found  by  multiplying 
by  the  span  length  of  each  panel: 


{L.}  -  p  u  {2S,r, } 
J  »  ®  j  j 


(A-12) 


Summing  the  lift  contributed  by  each  panel  gives  the  total  lift  for 
the  wing. 


*  -  l  ^  l  s  r 

j-i  J  j-i  J  J 


The  drag  forces  (D^  from  Fig  14)  are; 


(A-13) 


dj  ‘ 


sin  a  =  ijO 


(A-14) 
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<v  ■  Vi!  ■  p-“-  (2ViV 


(A-15) 


Figure  14.  Definition  of  Wing  Forces 


Resultant  moments  about  a  given  axis  can  be  calculated  from: 

N 

M  -  7  It.  (x.  -  x  )  (A-16) 

\  j-i  1  i  0 


jL  Li  (xj  -  V 


(A-17) 


This  method  produces  approximations  to  the  actual  aerodynamic 
forces.  The  accuracy  of  the  method  increases  as  the  number  of 
aerodynamic  panels  is  increased.  So  far,  the  circulation  singularity 
and  control  points,  (xj.yj)  and  (x^y^,  have  been  arbitrarily  located 
within  the  panels.  Dispersing  these  points  more  evenly  over  the  wing 
also  increases  the  accuracy  of  the  method.  Locating  the  circulation 
singularities  at  the  quarter-chord  mid-span  points,  and  the  control 
points  at  the  three  quarter-chord  mid-span  points  of  each  panel  has 
been  found  to  have  the  best  results. 

The  doublet  lattice  method  is  applicable  to  wings  of  arbitrary 
trapezoidal  shape.  The  aerodynamic  panels  need  not  be  rectangular, 
allowing  the  method  to  be  applied  to  tapered  and  swept  wings.  However, 
care  should  be  used  in  selecting  the  panel  size  and  the  locations  of 
'  >  the  circulation  singularities  and  control  points  for  swept  wings. 

The  aerodynamic  influence  coefficients  contain  a  possible  singularity 
if  the  x-coordinate  of  a  control  point  aligns  with  the  x-coordinate  of 
a  circulation  singularity  from  a  different  panel.  This  problem  can 
be  omitted  by  varying  the  span  of  each  box  with  respect  to  its  chord 
length.  Changing  either  the  number  of  span-wise  panels  or  the  number 
of  chord-wise  panels  will  usually  solve  the  problem. 

It  is  also  possible  to  vary  the  size  of  a  panel  with  respect  to 
other  panels.  This  allows  the  panelB  to  be  concentrated  close  to 
wing  tips  and  edges  where  the  circulation  may  change  drastically. 
Likewise,  panels  can  be  dispersed  over  the  central  portions  of  the 
wing  where  little  variation  is  expected. 

The  doublet  lattice  method  can  be  applied  to  wings  in  symmetric 
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flight  as  well  as  antisymmetric  flight.  Symmetry  can  be  accounted 
for  by  adding  the  effects  of  the  left  wing  to  the  aerodynamic  matrix. 
This  is  done  by  summing  the  effects  of  the  corresponding  aerodynamic 
panels  of  both  the  right  and  left  wings  on  an  element  of  the 
aerodynamic  matrix*  rather  than  the  effects  of  only  the  panel  on  the 
right  wing. 

The  doublet  lattice  method  is  easily  handled  with  a  computer  and 
also  lends  itself  to  other  purposes*  such  aB  the  eigenvalue  analysis 
developed  in  the  body  of  this  thesis. 


Appendix  B 

Listing  of  Computer  Programs 

This  appendix  has  two  parts.  The  first  part  is  a  listing  of  the 
NASTRAN  program  used  to  calculate  the  flexibility  coefficients.  The 
geometric  and  structural  data  listed  in  the  BULK  DATA  section  are  tor 
an  unawapt  wing  with  1002  of  the  fibers  oriented  along  the  sweep 
angle  (0  ■  0).  Varying  the  sweep  angle  affects  only  the  grid  point 
coordinates  and  varying  the  fiber  orientation  affects  only  the  last 
column  of  the  CQUAD2  element  cards. 

The  second  part  of  this  appendix  is  a  listing  of  the  Aerodynamic 
and  Eigenvalue  Analysis  program.  It  too,  is  for  an  unswept  wing.  The 
wing  can  be  "swept"  forward  by  changing  the  sweep  angle  on  the  DATA 
card.  The  wing  can  be  "rotated"  forward  by  also  changing  the  locations 
of  the  circulation  and  control  points  appropriately. 
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David  C.  Bannerroan,  born  4  May  1955,  lived  in  Frankenmuth, 

Michigan  for  the  first  18  years  of  his  life.  The  next  four  years 
were  spent  at  the  University  of  Michigan,  School  of  Engineering. 

He  graduated  with  a  Bachelor  of  Science  degree  in  Aerospace  Engineering 
in  April  of  1977.  The  same  day  he  received  his  commission  into  the 
United  States  Air  Force  through  the  ROTC  program.  David  then 
enrolled  in  the  School  of  Engineering,  Air  Force  Institute  of 
Technology  and  began  work  on  a  Master's  Degree  in  Aeronautical 
Engineering  in  October  of  the  same  year. 


UNCLASSIFIED 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1.  REPORT  MUMPER  a.  GOVT  ACCESSION  NO. 

AFIT/GAE/AA/7  8D- 1 

S.  RECIPIENT'!  CATALOG  NUMBER 

4.  TITLE  (tnd  Submit) 

A  STUDY  OF  THE  USE  OF  LAMINATED  COMPOSITE  WING 
STRUCTURES  FOR  CONTROLLING  THE  DIVERGENCE  SPEED 

OF  FORWARD  SWEPT  WINGS 

5.  TYPE  OF  REPORT  A  PERIOD  COVERED 

MS  Thesis 

4.  PERFORMING  OTO,  REPORT  NUMBER 

7.  AUTHORfaJ 

David  C.  Bannerman 

2nd  Lt  USAF 

t.  CONTRACT  OR  QRANT  NUMBERfiJ 

*.  PERFORMING  ORGANIZATION  NAME  ANO  ADDRESS 

Air  Force  Institute  of  Technology  (AFIT/ENY) 
Wright-Patterson  AFB,  Ohio  454.33 

10.  RRQORAM  ELEMtNT,  PROJECT,  T  ASK 
AREA  A  PORK  UNIT  NUMBERS 

11.  CONTROLLING  OFFICE  NAME  ANO  ADDRESS 

12.  REPORT  DATE 

December  1978 

IS.  NUMBER  OF  PAGE! 

72 

14.  MONlTORINO  AGENCY  name  «  AODRESSflf  dlllorent  from  Cotilrolllni  Olllct) 

1*.  SECURITY  CLASS,  (ot  Ihli  "port) 

UNCLASSIFIED 

"Tsi.  declassification 'downgrading 

SCHEDULE 

U.  DISTRIBUTION  STATEMENT  (o(  thl*  Report) 

Approved  for  public  release;  distribution  unlimited. 

17.  DISTRIBUTION  STATEMENT  (ol  tho  afulracf  anlarad  In  Block  30,  II  dlllottnl  Horn  Rapori; 

*UP,*LeMlNTA,,v  notes  Approved  for  public  release;  JoeeplKp.  lflp^s,  T^jor,  USAF 

IAW  AFR  190-17  Director  of  Information 

19.  KEY  WORDS  fConffriu*  on  rmvotto  //  n#c»ea«ry  tmd  Identify  by  block  numbor) 

Forward  Swept  Wings 

Composites 

Static  Stability 

Wing  Divergence 
v  , 

20.  Alfcjtf  RACT  (Conllnu*  on  rovotoo  oldo  It  necessary  ond  IdontUy  by  block  number) 

This  thesis  is  an  investigation  of  the  effects  of  composite  laminate 
orientation  on  the  divergence  problem  of  forward  swept  wings.  The  aerodynamic 
and  structural  properties  are  modeled  with  a  doublet  lattice  technique  and 
finite  element  method,  respectively.  The  static  equilibrium  equation  is 
applied  to  combine  the  two  models.  This  results  in  an  eigenvalue  problem  with 
the  divergence  dynamic  pressure  as  the  lowest  eigenvalue. 

A  computer  program  is  written  to  formulate  and  solve  the  eigenvalue  problem. 

It  is  applied  to  two  cases  of  example  wines.  For  the  first  case.  100%  of  the  , 

=*• 


DD  ,1 


FORM 
AN  71 


1473 


EDITION  OF  I  NOV  »»  IS  08S0LETE 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whan  Da  la  Entorod) 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  P  AOIfMian  AS  I  a  Smaratf; 


20(cont)  aw* — %FV 
fibers  are  varied  from  CP'  to  5u~ 
different  sweep  angles.  For  the 


5up  ahead-Qf  the  structural  swept  axis  for  five 
the  second  caseV-wlAt-J £7.  of  the  idbers  are 


varied  while  the  remaining  fibers  are  oriented  at  9Q  and  ffl  45^  to  the 
structural  sweep  axis.  /  ^ 

The  results  indicate  that  tj\&  divergence  problem  can  be  effectively 
controlled  and  possibly  eliminated  by  properly  orienting  the  composite 
laminates.  The  divergence  roeed  was  found  to  Increase  and  reach  a  maximum 
as  the  principal  fiber  direct ion-isy  swept  forward.  Maximum  values  occur  for 
fibers  oriented  between  15®''  and  ahead  of  the  structural  swept  axis  for 
wing  swept  forward  more  than 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  PAGEflWiFn  D Bnfud) 


